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Abstract 


-  r  This  work  follows  Capt,  Riddle 1 s  1989  thesis^which  compared 
conventional,  ,  and  controlled  guided  missiles  in  simulated  flyouts 
against  an  airborne  target.  The  simulation  results  did  not  clearly 
demonstrate  the  H„  and  H  controller  advantages.  In  this  studv  the  missile 

4  TP 

model  and  simulation  were  refined.  In  addition  to  simulation  runs,  , 

V  " 

performance,  noise  attenuation,  and  p  robustness  analyses  were  conducted. 

Results  showed  areas  where  improvement  was  needed  in  the  H„  and  H 

L  'T 

controllers  used  in  the  study.  An  Hn  and  H  controller  redesign  is 
recommended  using  the  p  synthesis  design  technique.  Additional  simulation 
runs  are  also  recommended  to  examine  whether  the  improved  controllers  reduce 

miss  distances.  /  "  "  . 

,  •  J  i  ,  '  '  '  * 
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ROBUSTNESS  OF  H„  AND  H  DESIGN  TECHNIQUES 
2  <r> 

IN  THE  DESIGN  OF  A  GUIDED  MISSILE  CONTROLLER 


I  Introduction 


1 . 1  Mot ivation 

The  objective  of  this  effort  is  to  show  the  advantages  of  H0  and  H 

2  cc 

controllers  through  their  application  in  a  surface-to-air  missile 
s i mu  1  at  ion . 

Capt  Riddle  began  this  project  with  his  1989  thesis,  "Design  of  a 

Robust.  Controller  for  an  Unstable  Nonminimum  Phase  Guided  Missile"  [1].  He 

constructed  a  missile  model  and  simulation,  formulated  a  design  technique 

for  H,.  and  H  controllers,  designed  an  H0,  H  ,  and  conventional  controller, 

and  measured  their  performance  with  miss  distance  results.  The  results 

showed  no  advantage  in  using  the  H„  and  H  controllers  over  the  conventional 

2  T 

cont  roi  ler . 

Building  upon  Capt,  Riddle’s  work,  the  approach  here  is  to  show  the 

advantages  of  H„  and  H  controllers  in  three  wavs:  1)  increase  simulation 
2  v 

fidelity,  2)  examine  limits  of  missile  performance,  and  3)  develop  a 
detailed  analysis  procedure  for  evaluating  the  controllers. 

This  study  shows  that  the  H£  controller  designed  by  Capt  Riddle  offers 
a  robustness  advantage,  but  could  be  improved  to  achieve  quicker  time 
response  and  noise  rejection  characteristics.  The  H^  controller  designed  by 
Capt  Riddle  was  found  to  be  a  faulty  design.  Unlike  the  intent  of  H^ 
optimization  to  provide  good  robustness  characteristics,  this  controller  was 
unable  to  stabilize  the  missile  when  high  frequency  unmodelled  dynamics  were 
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added.  The  missile  flyout  simulations  corresponded  well  with  the  controller 
analysis  predictions. 

The  analysis  and  simulation  developed  in  this  study  provide  a  very  good 
means  of  evaluating  advanced  controllers.  The  study  shows  that  the  and 
controllers  used  by  Capt  Riddle  were  not  optimally  designed. 

1.2  Overview 

To  assist  the  reader,  a  brief  mapping  of  the  document  contents  is 
presented  here.  Chapters  2,  3,  4  and  part  of  7  summarize  work  done  by  Capt 
Riddle.  Chapters  5  through  10  describe  the  progress  made  during  this  study. 

1.2.1  Missile  System  Description 

Chapter  2  describes  the  surface-to-air  missile  model  developed  by 
(apt  Riddle.  It  includes  the  high  frequency  unmodelled  dynamics  which  were 
added  in  this  study. 

1.2.2  Simulation  Description 

Chapter  3  describes  the  missile  flyout  simulation.  It  includes  a 

development  of  the  missile  simulation  transfer  function,  G  .  ,  and  the 

s  lm 

nonlinear  calculations  performed  by  the  simulation. 

1.2.3  Guidance  Law 

Since  the  guidance  law  plays  an  important  role  in  target  miss 
distances,  Chapter  4  is  devoted  to  a  description  of  the  guidance  command 
block  used  by  the  missile. 

1.2.4  Refining  the  Model  and  Simulation 

Three  areas  of  missile  model  and  simulation  improvement  are 
investigated  in  Chapter  5:  1)  reducing  simulation  time  step  size,  2)  adding 
high  frequency  unmodelled  dynamics,  and  3)  examining  plant  parameter 
var i at i ons . 
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1.2.5  Controller  Analysis 

It  is  difficult  to  measure  the  controllers  relative  capabilities  based 
on  target  miss  distances  so  a  detailed  controller  analysis  was  undertaken  in 
Chapter  6.  Performance,  noise  attenuation,  and  robustness  characteristics 
of  each  controller  are  examined. 

1.2.6  H„  and  H  Controller  Design 

2  CO 

The  controller  analysis  identifies  areas  where  improvements  can  be 

made  to  the  and  H  controllers.  Chapter  7  describes  the  process  used 

2  oo 

to  design  the  controllers  and  the  initial  values  for  weighting  filters 
chosen  for  controller  redesign. 

1.2.7  Simulation  Runs 

Chapter  8  describes  a  set  of  simulation  runs  which  were  made  with  the 
refinements  described  in  Chapter  5.  The  runs  are  limited  to  one  engagement 
scenario  to  reduce  scenario-dependent  miss  distance  variations. 

1.2.8  Conclusions  and  Recommendations 

Chapter  9  summarizes  what  was  learned  from  this  study  and  gives 
recommendations  for  follow-on  work. 

1.2.9  Missile  Transfer  Functions  To  State  Space 

An  important  part  of  the  missile  model  development  is  the  process  of 
converting  single-input  single-output  transfer  functions  into  state  space 
form  and  combining  them  to  create  a  state  space  model.  The  same  processes 
are  important  in  the  developemcnt  of  large  transfer  functions  in  state  space 
form  used  throughout  the  study.  To  aid  follow-on  efforts  and  to  confirm 
correctness  of  the  model  used,  the  developement  of  the  missile  model  is 
presented  in  Appendix  A. 
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1.2.10  Description  of  Code  Created  for  Study 

Again  with  the  intent  of  aiding  follow-on  efforts,  a  description  of  the 
features  of  the  computer  coding  developed  for  this  study  is  presented  in 
Appendix  B. 
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II  Missile  System  Description 


The  missile  used  in  this  study  is  the  same  generic  tail-control  missile 
used  by  Capt  Riddle  except  for  the  addition  high  frequency  unmodeled 
dynamics.  Details  of  the  origin  and  development  of  the  missile  model  can  be 
found  in  Capt  Riddle’s  thesis.  This  chapter  is  a  summary  of  the  model. 

Some  of  the  characteristics  of  the  nominal  missile  are  given  in  Table  2.1. 


Table  2.1  Some  Nominal  Missile  Characteristics 


Guidance 

Radar 

Cont ro 1 

Tail  Fin 

Speed 

Mach  3.5 

Mass 

127  kg 

Center  of  Gravity 

2.63  m  from  nose 

Inertia  (pitch  &  vaw) 

271.3  kg-m2 

Relationships  of  important  variables  are  shown  in  Figure  2.1. 

The  closed  loop  missile  system  can  be  represented  as  show.i  in  Figure 
2.2.  Its  components  will  be  described  in  detail  starting  with  the  nominal 
plant,  the  addition  of  the  high  frequency  unmodelled  dynamics,  the  guidance 
command  block,  and  finally  the  controller.  A  simplified  closed  loop  missile 
block  diagram  is  shown  in  Figure  2.3. 

2 . 1  Norn i nal  Plant 

The  assumption  is  made  that  the  missile  has  finished  boosting  and  is  in 
a  16  to  20  degree  climb.  Its  climb  rate  is  328  m/s  (1000  ft/s)  starting 
from  4921  meters  (15,000  ft)  altitude.  The  missile  model  is  two  dimensional 
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P  (x,y) 

M 


x  downrange 


Figure  2.1  Coordinate  Svsten 
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CLOSED  LOOP  MISSILE  SYSTEM 


c 


Figure  2.3  Simplified  Closed  Loop  Missile  Block  Diagram 


in  the  horizontal  plane;  however,  dynamic  pressure,  weight,  moment  of 
inertia  and  center  of  gravity  are  varied  as  if  the  missile  was  in  a 
sustained  propulsion  climb.  Dimensional  stability  derivatives  of  the 
missile  resulting  from  these  variations  are  given  by  the  following  equations 
and  shown  in  Figures  2.4  through  2.10.  These  equations  are  polynomials 
which  approximate  tabulated  values.  They  are  accurate  in  the  range  of  0  to 
8  seconds  of  missile  flyout  time. 


N  (t)  =  1.871  +  0. 202 1.  -  0.0262t2  +  0.0008t'3  (g’s/deg)  (2.1) 

(X 

N  (  t  )  =  1.016  +  0. 076t  -  0.0161t2  -  0.0004t.3  (g’s/deg)  (2.2) 

o 

M  ( t. )  =  0.0144  -  0. 3109t  -  O.Olllt2  +  0.0085t3  (deg/s2/deg)  (2.3) 

C? 

Mc(t)  =  -2.78  -  0 . 230t  +  0.044t2  -  O.OOlt3  (deg/s2/deg)  (2.4) 

0 

Vft)  =  1000  +  80 t  -  8t2  (m/s)  (2.5) 

M 

hu(t)  =  4921  +  328t  ( m )  (2.6) 

M 


The  nominal  flight  conditions  and  dimensional  stability  derivatives  used  in 
th<=>  controller  design  are  given  in  Table  2.2.  A  block  diagram  of  the 
resulting  nominal  plant  can  be  seen  in  Figure  2.11. 
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Table  2.2  Nominal  Missile  Parameters 


Altitude 

6096  m  (20,000  ft) 

Mach 

3.5 

VM 

1100  m/s 

N 

2.0  g’ s/deg 

« 

N6 

1.0  g’s/deg 

M 

-0.91  deg/s^/deg 

a 

M6 

-2.74  deg/s  /deg 

|RtI 

1500  m/s 

rt 

1000  n 

K 

3.3 

g 

NOMINAL  plant 


Figure  2.11  Missile  Model  Block  Diagram 
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The  resulting  state  space  representation  is  shown  in  Equations  (2.7)  through 


(2.14). 


x  =  A  x  +  B...U  +  B  w  +  B0C 

p  P  P  pH  P  pl2  p  p2  * 

v  =  C  x  +  D  , ,  u  +  D10w  +  D  „  r  +  n 

p  p  p  pi  1  p  pl2  p  p2  * 


where 


8 


X 

p 


8  (deg) 

0  (deg/sec) 

m 

0  (deg) 

D1 

cr  (deg) 

<P  (deg) 

0  (deg/sec) 
s 

0g  (deg) 

Xudl 

Xud2 

Xud3 


(2.7) 

(2.8) 


(2.9) 


u 

P 


6  (deg) 
com 

o  (deg/sec) 


(2.10) 


w 

P 


<*>T  (deg) 


(2.11) 
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y 


p 


An  (g’s) 

0  { deg/sec) 

ro 

0  (deg) 

0  (deg/sec) 
s 

6  (deg) 
s 

ud  (deg/sec  &  g’s) 


| 

(2.12) 


The  values  shown  for  £  and  q 


are  the  range  of  values  used  in  the  simulation. 


c  = 


+  (0  -  0.03) 


+  (0  -  3000)/Rt 
+  (0  -  0.5) 


turbulence  (deg/sec) 

target  glint  (deg/sec) 

power  supply  noise 
(deg/sec  &  deg) 


(2.13) 


n 


+  (0  -  1.0) 
+  (0  -  0.5) 
0 

+  (0  -  1.0) 
+  (0  -  0.5) 
0 


accelerometer  noise  (g’s) 

body  rate  gyro  noise  (deg/sec) 

no  measurement  noise 

rate  gyro  noise  (deg/sec) 

potentiometer  noise  (deg) 

nothing  added  to  high 
frequency  dynamics 


(2.14) 


The  state  space  values  for  the  plant  are  given  in  Equations  (2.15)  through 
(2.18).  Details  on  converting  the  missile  transfer  functions  to  state  space 
form  (’an  be  found  in  Appendix  A.  The  nominal  values  used  for  the  plant 
during  controller  design  can  be  seen  in  Equations  (2.19)  and  (2.20). 
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fin 


actuator 

0 

0 

0 

0 

0 

0 

0 

0 

0 

pole 

Mc 

0 

0 

-M 

0 

0 

0 

0 

0 

0 

6 

a 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

-L/Vu 

1 

0 

N  /V 

u  0 

0 

0 

0 

0 

0 

6  M 

a 

M 

-target 

target 

-target 

0 

0 

resolver 

0 

resolver  0 

resolver  0 

0 

0 

gain 

pole 

gain 

antenna 

0 

0 

0 

0 

0 

servo 

0 

0 

0 

0 

pole 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

high 

high 

high 

high 

freq  dyn 

0 

0 

0 

0 

0 

0 

freq 

freq 

f  req 

ga  i  u 

dyn 

dyn 

dyn 

pole 

pole 

pole 

-zero 

-zero 

h  i  gh 

high 

high 

freq  dyn 

0 

0 

0 

0 

0 

0 

0 

freq 

f  req 

ga  i  n 

dyn 

dyn 

pole 

pole 

-zero 

high 

high 

freq  dyn 

0 

0 

0 

0 

0 

0 

0 

0 

f  req 

gain 

pole 


(2.15) 
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c 

p 


BpU  Bpl2  Bp2 


10  0  2 

0  10  0 

0  0  0  0 

0  0  0  0 

0  0  0  0 

h  i  gh 

freq  dyn  0  0  0 

3a  i  n 


[  fin 

actuator  0 
gain 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

target 

1 

0 

0 

0 

0 

antenna 

*  resolver" 
gain 

0 

0 

1 

0 

servo 

gain 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

high 

high 

high 

freq 

f  req 

freq 

0 

0 

0 

dyn 

dyn 

dyn 

pole 

pole 

pole 

-zero 

-zero 

-zero 

o  o  •  0  ’  0  0  0 


D  .  ,  D  .  „  D  „ 
pi  1  pl2  p2 


0  0  *  0  *  0  0  0 

0  0  *  0  ‘  0  0  0 

0  0  *  0  *  0  0  0 

0  0  ■  0  ■  0  0  0 

0  0  •  0  •  0  0  0 


(2.16) 


(2.17) 


(2.18) 
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(2.19) 


G 

P 


G  = 

a 

< 

Bpll  Bpl 2  Bp2 

P 

LCP 

Bpll  Bpl 2  Bp2  _ 

-77 

0 

0 

0 

0 

0 

0 

0 

0 

0 

77 

0 

0 

0 

0 

1 

-2.74 

0 

0 

0.91 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-0.51 

1 

0 

-1.02 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

-20. 

7  0 

-5. 

26  0 

-20. 

7  0 

0 

0 

0 

0 

20. 

7  0 

0 

1 

0 

0 

0 

0 

0 

-104 

0 

0 

0 

0 

0 

94 

0 

0 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6.4 

0 

0 

0 

0 

0 

0 

-400 

-390 

-390 

0 

0 

0 

0 

0 

0 

6.4 

0 

0 

0 

0 

0 

0 

0 

-400 

-390 

0 

0 

0 

0 

0 

0 

6.4 

0 

0 

0 

0 

0 

0 

0 

0 

-400 

0 

0 

0 

0 

0 

0 

1 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

6.4 

0 

0 

0 

0 

0 

0 

-390 

-390 

-390 

0 

0 

0 

0 

0 

0 

_ 

J 

(2.20) 


2.2  Adding  The  High  Frequency  Unmodel  led  Dynamics 

Dif  unmodeled  high  frequency  dynamics  is  added  to  fl  and  A  by 

m  n 

multplving  v  bv  a  matrix  called  F  which  is  given  in  Equation  (2.21). 

'  p 


F 


10  0  0  0  1 
0  10  0  0  1 
0  0  1  0  0  0 
0  0  0  1  0  0 
0  0  0  0  1  0 


(2.21  ) 


2.3  Guidance  Command  Block 


The  guidance  command  block  can  be  seen  in  Figure  2.12.  Either  A  or 

ncom 

ft  can  be  used.  More  details  on  the  proportional  navigation  used  are  in 
mcom 
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GUIDE 


Chapter 
ca lied  T 


where 


The  561 . 


Figure  2.12  Guidance  Command  Block  Diagram 


5.  The  guidance  command  block  will  be  represented  by  a 
as  shown  in  Equations  (2.22)  through  (2.24). 


=  T 


y 


P 


v  =  A  or  0 
g  ncom  mcom 


1  0 

0  1 

T  =  0  0 

0  0 

0  0 


0  Kg]RT|/561 


0  0  0 

0  0  0 

1  0  0 

0  1  0 

0  0  1 

5  0  Kg|RT|/561.5  0 


5  in  Equation  (2.24)  would  be  replaced  by  1100  if  ^mcoin 


matrix 


(2.22) 


(2.23) 


(2.24) 


was  desired. 
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2.4  Controller 


Capt  Riddle’s  conventional  controller  can  be  seen  in  Figure  2.13. 


(Conv.  ) 
CONTROLLER 


Figure  2.13  Conventional  Controller  Block  Diagram 


This  controller  is  designed  for  use  with  negative  feedback.  Equations 
(2.25)  through  (2.2H)  show  how  the  controller  transfer  function  will  be 
represented . 

*  x  =  A  x  +  B  u  (2.25) 

c  c  c  c  c 

v  =  C  x  +  D  u  (2.26) 

'  c  c  c  c  c 

Dynamic  Controller 

(2.27) 

Static  Controller 

G  =  D  (2.28) 

c  c 
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Ill  Simulation  Description 


A  diagram  of  the  simulation  process  can  be  seen  in  Figure  3.1.  Each  block 
of  the  diagram  will  be  described  in  detail. 


*  Computation  shown  by  Equation  (3.8) 

Figure  3.1  Diagram  of  Simulation  Process 


Initial  conditions  used  in  the  missile  simulation  are  shown  in  Table 


3.1. 


3- 


Table  3.1  Simulation  Initial  Conditions 


Variable 

Value  Used  In 

Study 

PT  (  x ,  y ) 

0 

[  5000,  5000  ] 

meters 

VT 

470  m/s 

at 

-8.4  g’s 

-10  deg/s 

7? 

0 

0  deg 

PM  (  x  ,  y  ) 

Mo 

[  0,  0  ]  meters 

Vmo 

1000  m/s 

\ 

0  deg/s 

45  deg 

K 

0 

45  deg 

0 

0  deg/s 

rt 

0 

7071  meters 

\ 

-746  m/s 

\ 

0  deg 

0 

45  deg 

X 

P0 

0 

X 

co 

0 

^0 

0 

% 

0 

3.1  linearized  Scenario  System 

The  linearized  scenario  svstem,  G  used  in  the  missile  simulation  can 

sim 

be  seen  in  Figure  3.2. 
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(Simulation) 

LINEARIZED  SCENARIO  SYSTEM 


Two  wavs  to  obtain  v,  are  shown.  Both  yield  the  same  values.  Obtaining  v 

M  M 

from  A  was  used  in  the  simulation.  Capt  Riddle  used  (0  -  a)  to  obtain 

n  m  M 

in  his  simulation. 

The  closed  lo'-  ■  .ssile  system  state  space  representation  used  in  the 
missile  Simula1 .  is  shown  in  Equations  (3.1)  and  (3.2). 


Removing  u  and  u  through  substitution  using  u  =  -v  and  u  =  T  F  v 
pc  p  p  c  p 

results  in  Equations  (3.3)  and  (3.4). 
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Now  substituting  ^  =  R  y  where  R=[  .5100000]  and  integrating  to 
get  results  in  G  as  shown  in  Equations  (3.5)  through  (3.7). 
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G 

sim 


fA  -B  .  D  TFC  ] 
[  P  pll  c  pj 

-B  -  ,  C 
pll  c 

0 

(Bpi2-BpnDcTFDpi2) 

(Bp2-BPnDcTFy) 

(Bp2-BpllDcTF) 

B  TFC 

A 

0 

B  TFD  10 

B  TFD  „ 

B  TF 

c  p 

c 

c  pl2 

c  p2 

c 

RC 

P 

0 

o 

RDpl  2 

RDP2 

R 

0 

0 

1 

0 

0 

0 

(3.7) 


3.2  Time  Dependent  Target  Calculations 

Time  dependent  target  calculations  consist  of  the  following.  The 
target’s  turn  rate  (y^,)  is  calculated  for  a  given  target  speed  and  g’s  of 
acceleration.  For  a  constant  altitude,  turn  rate  is  given  by  Equation 
(3.8). 


f  AT  )  f 1801 

f9.8] 

0 

— 

1  A  1  -  1  7t  j 

V_ 

at^  -  1 

1  y  1  \  ) 

l  TJ 

k  > 

The  target  heading  is  calculated  at  each  time  step  through  Equation  (3.9). 


yT  -  tt  +  ip  (ts) 


(3.9) 


The  target’s  position  is  calculated  using  Equation 


PT(x,y)  = 


PT(x)  + 


vT(ts)cos(rT)),[pT(y)  t  VT(tJsin(,T)j 


(3. 10) 


3.3  Time  Dependent.  Missile  Calculations 

The  t  i  m^  dependent,  missile  calculations  consist  of  the  following.  The 


3-5 


plant  variations  and  missile  velocity  variation  which  were  described  in 

Chapter  2  are  updated  every  0.3  seconds  of  flyout  time  in  the  missile 

simulation.  The  heading  of  the  missile  (y„)  is  obtained  from  G  .  .  Missile 

1  sim 

position  is  calculated  each  time  step,  similar  to  target  position  with 
Equat ion  (3.11). 


PM(x,y)  = 


PM(X)  +  VM(ts)COSW)’fPM(y)  +  VM(tS)sin0 


(3.11) 


3.4  Nonlinear  Calculations  and  Time  Step  Differentiation 

Nonlinear  calculations  and  time  step  differentiation  consist  of  the 
following.  Angle  from  the  missile  to  the  target  with  respect  to  an  inertial 
reference  (<p^.)  is  calculated  from  missile  and  target  positions  with  Equation 
(3.12). 


<b 

'T 


(r 

arc  tan  j  p^  (  y ) 


-  PM<y 


)]/[PT(x)  ' 


1 1 

x’ii 


(3.12) 


* 

Range  from  missile  to  target  ( )  is  calculated  from  missile  and  target 
positions  with  Equation  (3.13). 


= 


[pTo 


PH,X|) 


-  PH(y>) 


(3.13) 


Range  rate  (R^,)  is  obtained  through  time  step  differentiation  (first  order 
backwards  difference  estimation)  as  shown  in  Equation  (3.14). 


R?  =  (rt ( i )  -  RT(i-l)]/ts 


(3.14) 
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3.5 


Noise  Generation 


The  derivation  of  the  equations  used  for  noise  generation  can  be  found 
in  [Ref  1].  The  r  and  r]  noise  generation  are  calculated  each  time  step  with 
Equations  (3.15)  through  (3.23). 

Turbu lence : 


=  .015  sinftj.)  +  sin(10  +  t^. )  +  .1  x  rand(l)j 


(3.15) 


Target  Glint: 


C.,  =  1500/R 


sin 


n 


t) 


rand (  1 

sin^lO  x  tj 


rand( 1 


3m) 


+  . 1  x  rand (  1  ) 


(3.16) 


Power  Supply  Noise: 


.25 


sin pOO  x  tf  ♦  -ranaT1 


3ttt) 


sin  ^2000  x  10  x  tf  +  — ^jyy)  +  .1  x  rand(l)j  (3.17) 


Accelerometer  Noise: 


=  .5 


s  in  j^40 


X  ^f  +  rand (  1  ) J 


sinj^400  x  tj.  +  ranj|(  |  ^  j  +  .1  x  rand(l)  (3.18) 


Body  Rate  Gyro  Noise: 

r>2  =  •25(siii(10  x  S-  +  + 

sin^lOO  x  tj,  +  rancHT)')  +  x  rand(l)j  (3.19) 
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No  Measurement  Noise: 


Hi  =  0 


(3.20) 


Antenna  Rate  Gyro  Noise: 


n4  -  .5  Sin(lO  »  tf  t 


sin^lOO  x  tf  +  -J— yj  +  .1  x  rand ( 1  )  (3.21) 


Antenna  Potentiometer  Noise: 


n. 

'.n 


=  .25 


i  i  n  f  t  „  + 


l  f  rand( 


tt) 


sinflO  x  t  +  - Jr— 

t  f  rand( 1  ) 


j  +  .  1  x  rand  (1)1 


1  x  rand ( 1 ) j  (3.22) 


Nothin?  Added  to  High  Frequency  Dynamics: 


n  ,  =  0 


3.23) 


The  rand(l)  function  generates  a  random  distribution  of  values  between  one 
and  minus  one.  A  typical  noise  generated  by  these  equations  during  a 
simulation  run  can  be  seen  in  Figures  3.3  and  3.4. 

3.6  End  Game 

Once  range  rate  is  no  longer  negative,  the  missile  is  going  away  from 
the  target,  so  end  game  calculations  are  made.  Calculation  of  the  distance 
of  closest  approach,  or  miss  distance,  is  as  follows.  Find  the  x  and  y 
components  of  the  range  from  missile  to  target  at  the  last  time  step  where 
Rt  i  0. 
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turbulence 


target  glint 


.  °'osh  AAA  a  a 

s  0 ! '  J  ^ V  v  ^  fi  A  A ! 

-o.o3.  v  V  y  v  V 


-0.04 


time,  sec 


time,  sec 


power  supply  noise 


0  l, 


-0.5  L 


time,  sec 


Figure  3.3  r  Noise  Inputs  Used  in  Simulate 
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(3.24) 


x1  =  PT(x)  -  PM(x) 
y,  =  PT(yl  -  PH(y) 


At  the  first  time  step  where  >  0,  define 


x  2  =  Pvp  (  x )  -  PM(x) 

y2  =  PT(y)  -  PM(y) 


(3.25) 


(3.26) 

(3.27) 


Write  the  \  and  y  range  components  in  terms  of  t  and  component  closure  rate, 
so  that 


x  =  Xj  +  t(dx/dt)  (3.28) 

y  =  >’j  +  t(dy/dt)  (3.29) 

Set  y  =  0  and  solve  for  t/dt. 

t/dt  =  -y^/dt  (3.30) 

Substitute  into  Equation  3.28  to  find  the  x  intercept  (x^ 

x.  =  x.  -  y.(dx/dy)  (3.31) 

int  1  1 


Using  the  same  process  find  v.  . 

int 


y .  .  =  v  -  x. ( dy/dx ) 
int  1  1 


(3.32) 
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The  relationship  between  x.  ,  y.  ,  and  miss  distance  (md)  is  shown  in 

int  i n  t 

Figure  3.5. 


Vsing  this  relationship  the  area  of  the  triangle  in  Figure  3.5  is  equal  to: 


Area  =  .5  x  x.  .  x  y.  , 
int  int 


Area  =  . 5  x  md  x 


2  ,  2 
x .  .  +  v . 

int  int 


(3.33) 


(3.34) 


Combine  Equations  (3.33)  and  (3.34)  to  solve  for  md. 


"d  =  (yint  x  xint)/ 


2  2 

x .  .  +  y . 

int  int 


(3.35) 


3-12 


IV  The  Guidance  Law 


Proportional  navigation  was  used  as  the  guidance  law.  Changing  the 
value  of  the  navigational  constant  in  the  guidance  law  plays  a  key  role  in 
miss  distance  results.  Breuer  (2:2-173]  gives  the  equation  for 
proportional  navigation  commands  as 


Ar  =  e  |Rt|  \t  (n/180)  m/s' 


(4.1) 


where  A  is  the  acceleration  of  the  missile  perpendicular  to  the  vector  from 
B 

the  missile  to  the  target  and  e  is  the  navigational  constant  used  by  Breuer. 

The  A  acceleration  will  maintain  a  constant  angle  between  the  target  and 
B 

the  missile  (A  )  bv  maintaining  A  at  zero.  The  value  chosen  for  r  is 
T  T 

related  to  how  tightly  the  missile  will  turn  to  achieve  a  new  desired 
head i ng . 

An  acceleration  normal  to  the  missile  velocity  vector  would  have  the 
following  component  of  acceleration  in  the  A  direction 


An  =  A  cos(AM)(9.8)  m/s 
Bn  M 


Substituting  into  Equation  (4.1)  and  solving  for  yields 


A  =  t /cost  A  )  |R_l  XT  (n/180  )/9. 8  g’i 
n  m  '  T 1  T 


(4.3) 


Fr<>m  Figure  2.1  it  can  be  seen  that 
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(4.4) 


X,j,  =  deg 

Therefore 


XT  =  <jf>T  -  yT  deg/s  (4.5) 

Substituting  this  result  into  Equation  (4.3)  yeilds 

=  t/cos(XM)  jR^I  (<f>T  -  yT)  (n/180)/9.8  g’s  (4.6) 

Tile  variables  of  Equation  (4.6)  must  now  be  transformed  into  quantities 
which  the  missile  can  measure.  Acceleration  normal  to  the  missile  velocity 
vector  can  be  approximated  by  accelerometer  readings  normal  to  the  missile 
body  axis  for  small  missile  angles  of  attack.  In  this  model  it  is  assumed 
that  the  accelerometer  readings  from  the  missile  are  the  same  as  the 
acceleration  normal  to  the  missile  velocity  vector,  or 

A  <s  A  (4.7) 

n  nm 


During  perfect  target  tracking 


a  +  e  =  <t>r 

m  s  r 


(4.8) 


so  if  the  seeker  is  tracking  well  can  be  approximated  with  8  and  8  , 

!  m  s 

which  are  measureable  with  rate  gyros. 
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(4.9) 


^m  +  6s  * 

The  terra  is  the  lead  angle  which  you  want  to  give  the  target.  This  angle 

will  depend  on  the  type  of  targets  you  are  engaging.  It  can  be  incorporated 

into  the  navigational  constant  called  K  . 

g 

K  =  c/cos (Xu)  (4.10) 

g  M 

Closure  rate,  | H ^ j  ,  is  measured  by  the  seeker  radar.  This  missile  is  unable 
to  measure  the  target  turn  rate  so  target  turn  rate  will  be  approximated  by 
0. 


iT  *  0 


(4.11  ) 


Substituting  Equations  (4.9)  through  (4.11)  into  Equation  (4.6)  yeilds  the 
form  of  the  guidance  commands  used  in  the  simulation. 


A  =  K  ! R_ I  (0  +  B  )  (n/180 )/9. 8  g’s 

ncom  g  1  T1  m  s 

=  K  j  R  I  (0  +  e  )  /561.5  g’s  (4.12) 

g  1  1  1  m  s 


This  ran  also  be  put  into  the  form  of  a  missile  heading  command  using  the 
rel at ionshi p 


9.8 


yu  =  (B  -  a)  deg/s 
M  m 


(4.13) 


Combining  Equations  (4.12)  and  (4.13)  using  the  nominal  value  for  of  1100 
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m/s  and  again  saying  angle  of  attack  of  the  missile  and  its  rate  are  small 
yields 


8  =  K  I R_ |  {8  +  6  )  /  1100  deg/s  (4.14) 

mcom  g  1  T 1  m  s 

This  is  the  guidance  command  used  by  Capt  Riddle. 

A  navigational  constant,  K  ,  equal  to  3.3  was  used  for  this  studv. 

g 

The  guidance  law  was  tested  by  running  the  simulation  with  no  missile 
dynamics;  that  is,  a  1100  m/s  point  moving  in  space  at  a  0.02  second 
simulation  time  step.  Target  miss  distances  were  only  a  few  centimeters 
caused  by  the  time  steping.  To  add  realism  the  simulation  limits  guidance 
gain  updates  to  once  every  0.3  second  and  commands  straight  ahead  flight 
near  end  game  when  R^/ 1 |  <=  0.6.  With  these  limitations  the  simulation 
with  no  missile  dynamics  still  came  within  2  meters  of  the  target  for  the 
engagement  chosen  for  this  study,  Figure  4.1.  This  gave  confidence  that  the 
miss  distance  which  resulted  would  be  due  almost  entirely  to  the  response 
lag  caused  by  the  missile  and  seeker  dynamics  and  not  due  to  the  limitations 
of  the  guidance  command  block. 

When  the  simulation  was  run  with  a  point  mass  responding  with  a  0.9 
second  lag  to  the  guidance  commands  miss  distances  corresponded  surprisingly 
well  with  those  resulting  from  the  simulation  run  with  full  missile  and 
seeker  dynamics,  Figure  4.2.  The  relationship  between  guidance  command 
execution  lag  and  resulting  miss  distance  for  the  chosen  engagement  can  be 
seen  in  Figure  4.3.  How  the  navigational  constant,  K  ,  effects  miss 

O 

distances  for  a  no  execution  lag  flyout  in  the  same  engagement  is  also  shown 
in  Figure  4.3. 
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Figure  4.1  Run  l,  Ideal  Missile  with  No  Execution  Lag 
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Figure  4.2  Run  2,  Ideal  Missile  with  0.9  Second  Execution  Lag 


ional  Constant  Effects 


When  running  the  simulations  with  full  missile  dynamics  it  was  found 
that  for  a  given  engagement  a  value  of  the  navigational  ~ould  be 

chosen  to  increase  the  missile  turn  rate  just  enough  to  offset  the  guidance 
command  execution  lag  at  the  point  of  end  game,  resulting  i.,  near  zero  miss 
distance.  That  was  equivalent  to  flying  a  timed  pattern  to  point  in  space 
where  the  target  would  be  when  the  missile  arrived.  It  is  entirely  scenario 
dependent  and  not  something  which  could  be  could  practically  applied  but  it 
does  point  out  that  using  miss  distance  for  a  measure  of  merit  in  judging 
competeing  controller  performance  can  sometimes  be  misleading.  That  is  why 
other  techniques  to  evaluate  controller  performance  were  used  along  with 
miss  distance. 
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V  Refining  the  Model  and  Simulation 


Three  areas  of  model  and  simulation  refinment  were  investigated. 
Simulations  were  run  with  a  smaller  time  step,  the  missile  parameter 
variations  were  checked,  and  high  frequency  unmodeled  dynamics  were  acHed  to 
the  missi le  model . 

Capt  Riddle  ran  his  simulations  with  a  0.02  second  time  step.  To 

increase  the  fidelity  of  the  simulation  a  selection  of  Capt  Riddle’s  runs 

were  made  with  a  0.02  and  a  0.001  second  time  step.  Runs  were  made  with  and 

without  noise  and  parameter  variations.  The  results  showed  there  was  no 

significant  difference  in  flyout  and  miss  distance  results.  Average  miss 

distances  w  -e  within  a  meter  of  each  other  for  the  two  time  steps. 

Capt  Riddle’s  simulation  runs  with  the  nominal  missile  performed  worse 

than  with  t  lie  missile  run  with  parameter  variations.  Though  this  seems 

unnatural  it  proved  to  be  resonable.  The  6  to  A  and  8  to  9  transfer 

com  n  com  m 

functions  are  shown  with  Eq  at  ions  (5.1)  and  (5.2)  and  in  Table  5.1.  The 
nominal  system,  which  is  roughly  the  system  at  4  seconds  of  flyout,  has  the 
nonminimum  phase  ze  o  and  the  unstable  pole  closest  together  making  it  the 
most  difficult  control  problem  to  solve.  Thus,  even  though  this  is  the 
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Transfer  Kunct  i  oiis 
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lit  roller  ties  i  i>'n  point,  it  is  not  unreasonable  t  iiat  variations  which  arc 
siot  to  control  per formed  tiottor.  During  this  investigation  it  was  found 
it  (apt  Riddle’s  Tuv  lor  series  approximations  for  -Nt/V'  and  -N  /V  were 

O  M  CX  PI 

t  very  ipiod  so  Taylor  series  approx imat  ions  were  replaced  with  the  actual 
vision  in  the  simulation  and  the  analysis. 


Hiilh  freipidiicy  uumodeled  dynamics  were  added  to  the  simulation  and  the 


missile  for  analysis  purposes.  They  are  not  included  for  controller  design. 


Through  their  addition  it  was  hoped  that  the  H„  and  H  controllers 

2  op 

robustness  advantages  could  be  seen.  Based  on  a  plot  of  measured  data, 
Reichert,  in  reference  [3:2371-2372],  adds  high  frequency  unmodeled  missile 
dynamics  to  A  using 


ud 

6“ 


100( s+10 ) 3 
( s+1000)3 


(5.3) 


This  was  adapted  for  this  study’s  missile  model.  The  high  frequency 

unmodeled  dvnamics  were  added  to  A  and  8  as 

n  m 

ud  6.4(s+10)3  ,, 

t~  =  - T"  (5 

°  (s+400) 

Equation  (5.4)  adequately  characterized  the  measured  data  and  the  smaller 
pole  value,  -400  versus  -1000,  makes  the  (sI-A)  3  matrix  inversion 
computationally  easier.  The  state  space  form  transfer  function  used  is 
shown  in  Equation  (5.5). 


-400 

0 

0 

-390 

-390 

-400 

0 

-390 

-390 

-390 

-400 

-390 

6.4 

6.4 

6.4 

6.4 

Appendix  A  shows  how  this  was  integrated  into  the  plant  state  space 
matrices.  Chapter  8  shows  the  effects  the  high  frequency  unmodeled 
dynamics  had  on  the  simulations. 
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VI  Controller  Analysis 


Because  Capt  Riddle’s  simulation  miss  distance  results  did  not  clearly 
differentiate  desirability  of  the  three  controllers,  other  analysis 
techniques  to  supplement  miss  distances  were  deemed  necessary.  The 
controllers  were  evaluated  by  evaluating  the  closed  loop  missile  in  three 
areas:  performance,  noise  attenuation,  and  robustness. 

6.1  Performance  Analysis 

The  missile’s  ability  to  carry  out  commanded  inputs  with  minimal  lag 
was  evaluated  by  taking  the  guidance  command  block  derived  in  Chapter  6 


A  =  K  |RJ  /  561.5  g’s  (6.1) 

ncom  g  1  I  1  1 

converting  A  to  turn  rate 

ncom 

|  Ry  1  <i>T  /  1100  deg/s  (6.2) 

and  input ing  a  step  command  for  ^  using  the  nominal  values  for  | R^ | ,  and 
K  shown  in  Table  6.1. 

o 

Table  6.1  Nominal  Values  Used  in  Performance  Evaluation 
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The  value  of  10  deg/s  for  0^  was  chosen  after  running  some  simulated 

engagements  to  find  the  upper  limit  which  0^  reached  outside  of  the  endgame. 

The  result  is  a  commanded  acceleration  normal  to  the  missile  velocity  vector 

of  88  g’s  or  a  turn  rate  of  45  deg/s.  The  0^  to  transfer  function  was 

found  for  the  nominal  missile  system.  The  yu  time  response  was  examined  for 

M 

a  0^.  step  input  of  10  deg/s.  This  was  compared  to  the  commanded  turn  rate 

of  45  deg/s  as  a  measure  of  controller  performance.  Singular  value  plots 

for  this  transfer  function  were  also  examined. 

To  find  the  0  to  y  transfer  function,  first  calculate  the  closed  loop 

transfer  function  for  the  nominal  missile,  G  .  Figure  6.1  shows  the  block 

m 

diagram  and  Equation  (6.3)  shows  the  transfer  function  in  state  space  form 

of  G  . 
m 


Figure  6.1  Simplified  Closed  Loop  Missile  Block  Diagram 


G  = 
m 

(A  -B  ,.D  TFC  ) 
p  pile  p 

B  TFC 
c  p 

-B  ,.C 
pi  1  c 

A 

c 

Bpl  2 

0 

C 

P 

0 

0 

Next ,  an  integrator  is  put  in  front  of  the  input  and  R  =  [.51  00000) 
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is  multiplied  by  the  output  to  obtain  the  transfer  function  for  <j>  to  y  , 

T  M 

G  .  The  block  diagram  for  G  is  shown  in  Figure  6.2  and  its  state  space 
s  s 

transfer  function  is  shown  in  Equation  (6.4). 


Figure  6.2  Gg  System  Block  Diagram 
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The  results  of  this  analysis  is  shown  in  Figures  6.3  through  6.6.  It 
is  done  with  and  without  the  unmodeled  dynamics  added. 

The  closed  loop  missile  time  responses  with  the  the  and 
controllers  are  better  damped  than  the  conventional  controller  but  without  a 
reduction  in  lag  time.  The  only  difference  in  time  responses  when  the  high 

frequency  unmodeled  dvnamics  are  added  is  that  the  H  controlled  missile 

oc 

goes  unstable.  This  is  contrary  to  optimization  theory  indicating  that 
this  is  an  improperly  formulated  H  controller. 

'  'T> 

The  closed  loop  missile  singular  value  plots  should  have  a  flat  gain  of 
4.5  (13  (IB)  at  low  frequencies  with  a  gain  drop  off  at  high  frequencies. 

The  singular  value  plots  of  t he  closed  loop  missile  with  the  and 
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Compensator: 

Riddle  Conv  — — 

Riddle  H-2  - 

Riddle  H-inf  .... 

- 1 1  J.  .1.. 1.1.1  i - 1 1 1 L  11.1 - 1 1 1—1  .l.i.  1.1.1 - J 1 I L-I  1  1  1  J _ I t l — L  .  L  J  111 1 l L 

IQ'2  10'1  10°  10*  102 


Frequency  (rad/sec) 


Figure  6. 4  Singular  Value  Plots  without  Unmodelled  Dynamics 


Singular  Values  dB 


Singular  Value  Plots  of  System,  RTdot  =  1500,  undyn 


Figure  6.6  Gg  Singular  Value  Plots  vith  Unmodelled  Dynamics 
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controllers  are  better  than  those  with  the  conventional  controller  at  high 

frequencies  but  are  lacking  in  gain  at  the  lower  frequencies.  This 

indicates  the  need  for  some  filter  shifting  during  the  H0  and  H  controller 

i.  <s> 

design  process.  Adding  high  frequency  unmodelled  dynamics  only  slightly 
effected  the  singular  value  plots  at  high  frequencies. 

6.2  Noise  Attenuation  Analysis 

For  this  analysis  the  same  nominal  system  is  used  as  in  the  performance 
analysis  except  this  time  the  q  to  y^  and  £  to  y^  transfer  functions  are 
found.  Singular  value  plots  of  the  transfer  functions  are  analyzed  to  see 
how  well  noise  inputs  are  attenuated.  Again  cases  with  and  without  high 
frequency  unmodeled  dvnmics  are  examined. 

The  block  diagram  for  the  rj  to  y  transfer  function  is  shown  in  Figure 

6.7. 


Figure  6.7 


q  to  y^  Block  Diagram, 


G 

n 


N  is  a  matrix  which  removes  q^  and  q^  inputs  since  there  are  none  in  the 
simulation.  N  is  given  in  Equation  (6.5). 


N 


10  0  0 
0  10  0 
0  0  0  0 
0  0  10 
0  0  0  1 
0  0  0  0 


(6.5) 
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The  state  space  transfer  function  for  r  to  is  shown  in  Equation  (6.6). 


(A  -B  ,,D  TFC  ) 
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RC 

L  P 

0 

0 

The  block  diagram  for  the  £  to  y^  transfer  function  is  shown  in  Figure 

6.8. 


Figure  6.8 


£  to  Block  Diagram, 


G 


C 


The  state  space  transfer  function  for  £  to  y^  is  shown  in  Equation  (6.7). 


A  0 

P 

0 

BP2 

Cc  = 

-B  ,.D  TFC  (A  -B  .  .D  TFC  ) 

pile  p  p  pile  p 

-B  ..C 
pi  1  c 

-BpllDcTFDp2 

B  TFC  B  TFC 

c  p  c  p 

A 

c 

B  TFD  .. 
c  p2 

RC  RC 

P  P 

0 

RDp2 

(6.7) 


The  noise  analysis  results  are  shown  in  Figures  6.9  through  6.12.  Again  13 
dB  would  be  zero  attenuation.  In  general  both  the  r]  and  £  attenuation  are 
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Magnitude  Jb 


10-3  10-2  1 0-*  10°  101  10s  103 

Frequency  (rad/sec) 


Figure  6.9 


Singular 


Value  Plots  without  Unmodelled  Dynamics 
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System  Noise  Analysis,  Zeta  Inputs,  RTdot  =  1500,  undyn 


Frequency  (rad/sec) 


ure  6.11  G  Singular  Value  Plots  with  Unmodelled  Dynamics 
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P 


System  Noise  Analysis,  Eta  Inputs,  RTdot  =  1500,  undyn 


Figure  6.12  G  Singular  Value  Plots  with  Unmodelled  Dvnamics 

n 


6-13 


better  with  the  conventional  controller  than  with  the  H0  and  H  controllers 

2 

indicating  a  need  to  change  the  filters  used  during  and  design. 

6.3  p  Robustness  Analysis 

The  jj  robustness  analysis  was  performed  as  described  in  reference  [4]. 
The  code  from  this  reference  is  also  used.  The  first  step  in  the  p  analysis 
is  to  construct  an  M  matrix.  The  form  of  the  M  matrix  is  shown  in  Figures 
6.13  and  6.14. 


Figure  6.13  Formation  of  M  Transfer  Function  for  p-Analysis 
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Figure  6.14  M  Transfer  Function  for  p-Analysis 
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Then  a  scale  value  is  chosen  to  indicate  the  extent  of  the  variation.  Plant 
variations  considered  are  those  due  to  aerodynamic  changes  with  altitude  and 
missile  weight  changes  as  described  in  Chapter  2.  Variations  in  the 
guidance  law  gain  are  also  considered.  The  guidance  law  is  allowed  to  vary 
by  changing  |  R,j,  |  from  0  to  2000  m/s.  Again  a  nominal  value  at  the  midpoint 
of  the  range  is  chosen  with  scaling  inputs  indicating  the  extent  of  the 
variation.  The  scaling  information  is  contained  in  a  matrix  called  PSC.4L 
shown  in  Table  6.2. 


Table  6.2  Plant  Scaling  Matrix  Used  in  p  Analysis 
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All  the  contro1  lers  maintain  a  value  of  u  ,  below  one  for  the  case 

'real 

without  high  frequency  unmodelled  dynamics,  indicating  the  closed-loop 
systems  are  stable.  When  high  frequency  unmodelled  dynamics  are  added,  only 
the  system  with  the  controller  remains  stable  for  high  frequency  inputs. 
The  performance  of  the  H  controller  is  worse  than  expected  again  indicating 
an  error  was  made  during  the  formulation  of  the  controller  design. 
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VII  H,.  and  H  Controller  Design 

Z  'T' 

Due  to  lack  of  time  an  improved  controller  design  was  not  completed. 
However,  a  brief  background  of  the  design  process  used  will  be  described  in 
this  chapter. 

The  appeal  of  and  controller  designs  is  that  they  minimize 
disturbance  effects,  w,  on  outputs,  z.  This  is  shown  in  block  diagram  form 
in  Fi gnre  7,1. 


Figure  7.1 


Plant  Block  Diagram  for  H,.  and  H  Construction 

2  a 


The  variable  definitions  as  applied  to  the  problem  in  this  study  are  given 
l n  Tab  1 e  7.1. 

The  H  design  objective  is  to  minimize  the  2-norm  shown  in  Equation 
(7.1). 
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The  2-norm  is  related  to  the  area  under  a  singular  value  plot. 
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i p 1 e- i nput/mu 1 1 i pi  e-out  put  designs  weighting  filters  are  used 


to  tailor  the  relative  importance  and  frequency  ranges  of  different  inputs 
and  outputs  in  the  controller  design.  Input  weighting  filters  are  chosen  to 
look  like  expected  disturbance  or  command  inputs.  Output  weighting  filters 
or  penalty  functions  are  chosen  with  a  frequency  band  to  match  the  desired 
output.  The  higher  the  magnitude  of  a  penalty  function  for  a  commanded 
input  the  harder  the  controller  will  try  to  obtain  that  commanded  input 
relative  to  the  other  inputs. 

The  block  diagram  form  of  the  transfer  function  G,  shown  in  Figure  7.1, 
for  this  problem  with  the  required  weighting  filters  is  shown  in  Figure  7.2. 


Figure  7.2 

Block  Diagram  Showing  Weighting  Filters  Used  In 
Construction  of  H9  and  H  Controllers 


The  compensator  in  Figure  7,2  is  considered  open.  The  weighting  filters  are 


shown  in  Equations  (7.3)  through  (’ 

1.7). 
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The  state  space  form  of  G  is  shown  in  Equations  (7.8)  through  (7.16). 
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The  derivation  of  these  matrices  is  similar  to  the  derivation  of  the  missile 

state  space  transfer  function  found  in  Appendix  A. 

The  Bode  magnitude  plots  of  the  weighting  filters  chosen  are  shown  in 

Figures  7.3  through  7.7.  The  magnitudes  of  these  filters  were  then  scaled 

to  produce  the  desired  affect.  Once  the  filters  were  chosen,  automated 

and  design  routines  were  employed  to  get  the  controllers.  The 

information  to  write  the  Hn  and  H  design  routines  can  be  found  in  [5]. 

2  -v 

Finally,  the  controllers  were  multiplied  by  a  -H  matrix  to  convert  v 

g 

outputs  of  the  guidance  block  into  the  ( r-v )  inputs  expected  by  the 
controller  and  to  convert  the  controller  from  one  designed  for  positive 
feedback  to  one  that  uses  negative  feedback,  as  shown  in  Figure  7.8. 

The  values  chosen  for  H  are  shown  in  Equation  (7.17). 
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The  resulting  compensator  is  run  through  the  analysis  described  in 
Chapter  6  and  then  the  scaling  on  the  filters  is  adjusted  to  get  the  desired 
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Figure  7.8  FL  and  H  Controller  Adaptation 

2  T> 


improvements.  The  process  is  difficult  because  there  are  many  filters 
involved  and  a  change  in  one  can  often  have  unexpected  effects  on  the 
others.  A  working  controller  is  achieved. 
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VIII  Simulation  Runs 


One  scenario  is  selected  for  the  simulation  runs  to  eliminate  scenario 
dependent  variations  in  miss  distance  and  bring  out  variations  due  to  *he 
missile  controller  performance.  A  turning  tail  chase  scenario  is  selected 
with  the  target  pulling  8.4  g’s  at  Mach  1.5.  It  was  hoped  that  the  added 
missile  flyout  time  of  a  tail  chase  would  better  differentiate  controller 
performance.  The  8.4  g  target  turn  is  chosen  because  the  navigational 
constant  { >  of  3.3  used  in  this  study  and  Capt  Riddle’s  is  best  suited  for 
highly  maneuverable  targets.  The  run  matrix  chosen  is  shown  in  Table  8.1. 

8.1  1 dea 1  F 1 yout 

First  the  simulation  is  run  with  no  missile  dynamics  to  confirm  the 
ef f feet iveness  of  the  guidance  command  block,  Run  1,  Figure  5.1.  Next,  the 
simulation  is  run  again  with  no  missile  dynamics  but  with  a  0.9  second  lag 
in  carrying  out  the  guidance  commands  as  a  point  of  comparison  for  the  lags 
caused  by  the  actual  missile  dynamics,  Run  2,  Figure  4.2. 

8.2  Baseline  Runs 

A  baseline  run  is  made  for  each  controller  with  the  nominal  plant,  no 
h  i  sr  ti  frequency  unraodelled  missile  dynamics,  and  no  disturbances  (noise),  Run 
3  through  Run  5. 

8.3  Plant  Variation  Effects 

The  effects  of  plant  variations  alone  are  examined  for  each  controller 
with  Runs  6  through  8.  Again  the  miss  distances  shown  in  Table  8.1  are 
improved  with  the  plant  variations  compared  to  the  nominal  plant  as  was 
explained  in  Chapter  5. 
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Table  8.1  Simulation  Run  Selection 


Run 

Number 


C’ont  roller 

Features 

Time  Step 
( seconds ) 

Miss  Distance 
( meters ) 

guide 

no  lag 

0.02 

1.5G 

gu  ide 

0.9  sec  lag 

35.04 

conv . 

Basel ine 

49.00 

H2 

pnom  nodvn  nodis 

0.02 

47.35 

H 

Y 

102.30 

conv . 

F’lant  Variations 

40.90 

H2 

pvar  nodvn  nodis 

0.02 

36.67 

H 

72.49 

conv . 

Unmodel  led 

48.81 

Dynam i cs 

H2 

pnom  undyn  nodis 

0.001 

47.56 

H 

. 

185.50 

conv . 

Noi  se 

48.89,  48.89 
48.89 

H2 

pnom  nodvn  noise 

0.001 

46.84,  47.26 
46.84 

H 

\ 

110.40,111.40 

.. 

110.40 

conv  . 

Combined  Effects 

40.85,  40.61 

H  V 

pvar  undyn  noise 

0.001 

37.45,  35.79 

H 

4 

48.02.646.30 

Hi'jh  Frequency  Unmodelled  Dynamics  Effects 


Th>  ■  f  feet  s  of  high  frequency  unmodelled  dynamics  alone  were  examined 
for  each  controller  in  Runs  9  through  11  The  step  size  of  0.001  second  was 
used  to  make  sure  the  simulation  captured  these  effects.  As  predicted  by 
the  ;i-ana 1\ sis  in  Chapter  6,  the  missile  with  the  "psuedo  "  controller  can 
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hardlv  he  called  stable,  Figure  8.1.  The  p-analysis  also  predicted  the 
missile  with  the  conventional  controller  would  be  unstable  with  unmodelled 
dynamics.  The  commanded  inputs  achieved  in  this  engagement  may  not  have 
been  high  enough  or  of  the  right  frequency  to  bring  out  the  predicted 
instability.  Each  time  the  fin  of  the  missile  is  jerked  quickly  the  missile 
rings  from  the  high  frequency  unmodelled  dynamics  shown  in  Figure  8.2  taken 
from  Run  9  and  Figure  8.3  taken  from  Run  10. 

8.5  Noise  Effects 

The  effects  of  noise  added  to  the  baseline  runs  are  examined  with  Runs 
12  through  20.  The  noise  simulations  are  run  more  than  once  to  see  how  much 
the  simulation  results  change  each  time  due  to  the  randomness  of  the  noise 
inputs.  Agian  a  0.001  second  time  step  is  chosen  to  ensure  high  frequency 
noise  effects  can  be  seen. 

The  miss  distances  with  the  conventional  and  controller  change  vary 
little  from  their  baseline  values  with  added  noise.  Miss  distance  values 
with  the  H  controller  increase  somewhat  over  its  baseline  values.  The 

A 

noise  attenuation  ability  of  the  controllers  is  shown  in  Figures  8.4 
through  8.8.  The  missile  with  the  H?  controller  has  the  least  noise  in 
O' <mmand*-d  turn  rate  tracking  followed  closely  by  the  missile  with  the 
conventional  controller.  The  missile  with  the  H  controller  has  about  three 


times  the  variations  in  commanded  turn  rate  as  the  other  two.  Actual  turn 

rate  variations  are  about  equal  for  the  three  controllers  in  high 

frequencies  but  in  low  frequency  deviations,  the  conventional  controller 

clearlv  does  better  followed  by  the  H„  and  the  H  controllers. 

2  t- 

8.8  Combined  Effects 

Finally  the  combined  effects  of  plant  variations,  high  frequency 
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Figure  8.1  Run  11,  Unmodelled  Dynamics  Effects  on 
Missile  with  H  Controller 


High  Frequency  Unmodelled  Dynamics 


8-5 


Figure  8.2  Run  9,  Generating  High  Frequency  Unmodelled  Dynamics, 
Conventional  Controller,  Nominal  Missile 


High  Frequency  IJnmodelled  Dynamics 


High  Frequency  Unmodelled  Dynamics, 
,  Nominal  Missile 


Desired,  Commanded,  and  Actual!  Missile  Turn  Rate 


Figure  8,4  Run  14,  Turn  Rate  Responses, 
Conventional  Controller,  Nominal  Missile  with  Added  Noise 


Desired,  Commanded,  and  Actual  Missile  Turn  Rate 


F'it»n re  H.5  Run  17,  Turn  Rate  Responses, 
Controller,  Nominal  Missile  with  Added  Noise 


Desired,  Commanded,  and  Actual  Missile  Turn  Rate 
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Ki ?ure  8.8  Run  20,  Turn  Kate  Responses, 
Controller,  Nominal  Missile  with  Added  Noise 


unmodelled  dynamics,  and  noise  inputs  are  examined  for  each  simulation  with 
Runs  21  through  26.  The  miss  distances  with  the  conventional  and 
controllers  are  similar  to  their  baseline  runs  with  plant  variations.  Again 
the  missile  with  the  controller  goes  unstable  though  it  still  manages  to 
get  near  the  target  in  Run  25  but  not  Run  26. 

Comparing  Figures  8.7  ana  8.8  it  can  be  seen  that  commanded  turn  rate 
has  less  variations  with  the  Hg  controller  than  with  the  conventional 
controller.  It  can  also  be  seen  that  the  conventional  controller  has  less 
low  frequency  turn  rate  deviations  and  slightly  more  high  frequecv 
deviations  than  with  the  controller.  Both  match  the  0.9  second  lag 
approximation  fairly  well. 

Adding  noise  brought  out  more  of  the  high  frequency  unmodelled  dynmics 
with  all  the  controllers.  This  effect  can  be  seen  in  Figure  8.9  with  the 
conventional  controller  and  Figure  8.10  with  the  controller. 
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Figure  8.7  Run  22,  Turn  Rate  Responses, 
Conventional  Controller,  Combined  Effects 
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Figure  8.8  Run  24,  Turn  Rate  Responses, 
H„  Controller,  Combined  Effects 


Frequency  Unmodelled  Dynamics 
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re  H.9  Run  21,  Generating  High  Frequency  Unmodelled  Dynamics 
Conventional  Controller,  Combined  Effects 


Frequency  Unmodelled  Dynamics 


Generating  High  Frequency  Unmodelled  Dynamics, 
Controller,  Combined  Effects 


IX  Conclusions  And  Recommendations 


9. 1  Conclusions 

Care  must  be  taken  when  judging  missile  controllers  based  on  target 
miss  distances  obtained  from  missile  flyout  simulations  even  though  low  miss 
distance  is  the  desired  real  world  end  result.  This  is  due  to  the  fact  that 
the  chosen  engagement  scenario  and  navigational  constant  also  drive  miss 
distance  values.  Reducing  command  execution  lag  time  caused  by  the  missile 
dynamics  has  the  greatest  impact  on  reducing  miss  distances.  This  lag  is 
the  combination  of  the  seeker  tracking  lag  and  lag  in  pointing  the  missile 
velocity  vector.  A  plot  of  desired,  commanded,  actual,  and  desired  with  lag 
missile  turn  rates  obtained  during  the  missile  flyout  simulation  proved  very 
useful  in  evaluating  controller  performance.  The  controller  can  compensate 
for  the  missile  dynamics,  but  there  is  a  physical  limit  to  what  the  missile 
seeker  combination  can  do  as  is  well  shown  by  Nesline,  William  and  Nesline 
in  reference  [ 6  J . 

The  results  of  the  controller  analyses  correlated  well  with  the 
simulation  runs.  The  controller  did  not  perform  as  it  should  have  and 
must  have  been  incorrectly  formulated  by  Capt  Riddle.  The  controller 
passed  the  robustness  tests,  but  other  analyses  showed  its  noise  attenuation 
and  time  response  could  be  improved.  A  redesign  of  these  controllers  was 
initiated,  but  time  limitations  prohibited  the  completion  of  an  improved 
controller  design.  Correlating  the  choice  of  weighting  filters  in  the 
controller  design  with  a  given  controller  performance  result  was  more 
difficult  than  anticipated. 
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9.2  Recommendations 


This  study  was  a  follow  on  to  Capt  Riddle’s  1989  thesis  work.  The 
model  and  simulation  used  is  therefore  in  its  second  generation.  Because  of 
this,  the  documentation  is  much  improved  and  the  computer  code,  despite  its 
added  features,  is  more  readable.  There  are  still  more  benefits  to  be 
gained  by  continuing  this  study.  It  is  suggested  that  the  H0  and  H 

Cj  00 

controllers  be  redesigned  based  on  the  information  gained  from  the 
controller  analyses.  An  iterative  design  technique  using  p-analysis  could 
ensure  the  full  potential  of  the  controllers  would  be  reached.  Simulation 
runs  would  then  determine  whether  or  not  reduced  miss  distances  could  be 
acheived  with  this  missile,  p-synthesis  is  another  option  for  designing  the 
controllers  17]. 

Another  suggested  area  for  follow-on  work  is  to  improve  the  controller 
analysis  by  using  p-analysis  for  noise  attenuation  and  performance  as  well 
as  robustness.  The  performance,  noise  attenuation,  and  robustness  analyses 
could  then  be  combined  into  one  p— analysis.  The  p-analvsis  computer  code 
(41  used  in  this  study  was  unable  to  handle  performance  and  noise  analysis, 
which  is  why  time  responses  and  singular  value  plots  of  nominal  transfer 
functions  were  used  instead.  The  p  robustnes,  performance,  and  noise 
anaylsis  block  diagrams  are  shown  in  Figures  9.1  through  9.3  respectively. 
The  combined  performance,  noise  attenuation,  and  robustness  p  analysis  is 
shown  in  Figure  9.4. 
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Figure  9.4 


y  Combined  Performance,  Noise,  And  Robustness  Analysis 
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Appendix  A:  Missile  Transfer  Functions  to  State  Space 


Converting  the  missile  transfer  functions  to  state  space  in  Chapter  2, 
forming  the  transfer  function  for  the  missile  simulation  in  Chapter  3,  and 
forming  the  transfer  functions  for  the  controller  analysis  in  Chapter  6  all 
make  use  of  the  following  equations.  "Transfer  function  to  state  space 
operation"  is  shown  with  Equations  ( A . 1 )  through  (A.4). 


-den( 2: N ) 

A  =  ( A .  1  ) 

I(N-2,N-1) 


where 


B  =  KN-1,1) 

C  =  num(:,2:N)  -  nun(:,l)  den(2:N) 
D  =  num( : ,  1 ) 


matrix  of  the  numerator  coefficients  of  the 
transfer  function 

vector  of  the  transfer  fuction’s  denominator 
coefficients 


I ( m,n ) 


identity  matrix  with  dimensions  m  rows  by  n 
columns 


power  of  transfer  function’s  characteristic 
equation  (denominator)  plus  one 


:  =  through  or  all 


(A. 2) 
(A. 3) 
(A.4) 


Transfer  function  multiplication  or  series  operation  is  shown  in  Figure  A.l 
and  Equations  (A. 5)  through  (A.l).  These  equations  may  not  yield  a  transfer 
function  with  minimal  states. 
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Transfer  function  addition  or  parallel  operation  is  shown  in  Equations  (A. 8) 
through  (A. 9)  and  in  Figure  A. 2.  Again  a  transfer  function  with  minimal 
states  may  not  be  produced  with  these  equations. 


C  =  G1  +  C2 


(A. 8) 


A-2 


A1 

0 

B1 

G  = 

0 

A2 

B2 

,C1 

C2 

VD2 

H  G 


Figure  A. 2  Parallel  Operation 


Since  the  transfer  functions  used  all  started  as  single  input  single  output 
eiiuat ions,  the  transpose  relationship  in  Equation  (A. 10)  can  also  be  used. 
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Equations  (A.  11)  through  (A. 17)  briefly  go  through  the  formation  of  the 
missile  state  space  transfer  function  formation  from  Figure  2.11  to  Equation 
(2.20). 
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Equations  (A. 11)  through  (A. 17)  are  combined  to  form  the  plant  state  space 
transfer  function  seen  in  Equation  (2.20). 
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Appendix  B:  Description  of  Code  Created  For  Study 


Pro-Mat  lab,  a  control  design  program,  was  used  for  this  study. 
Programming  can  be  done  within  Pro-Matlab  using  m  files.  Capt  Riddle’s  code 
was  used  as  a  starting  point  since  it  performed  some  of  the  same  functions 
needed  for  this  study.  The  code  was  entirely  rewritten  in  an  attempt  to 
make  it  easy  to  understand,  to  identify  errors,  and  to  modify  it.  It  is 
heavily  commented.  Much  of  this  document  was  taken  directly  from  the 
commented  code.  The  computer  code  consists  of  five  "m"  files,  four  used  in 
controller  design  and  evaluation  shown  in  Figure  B.l  and  the  fifth  used  to 
study  the  missile’s  navigational  command  block.  Though  the  code  is  readable 
and  self  explanitory,  some  of  its  features  will  be  described  here. 

Missile. m  is  a  function  which  contains  all  the  missile  state  space 
information  described  in  Chapter  2.  It  also  has  the  equations  needed  to 
compute  the  transfer  functions  used  in  the  controller  analysis.  Mrform.m  is 
called  by  missile. m  to  form  the  M  matrix  for  analysis. 

The  weighting  filters  are  chosen  by  the  designer  in  build. m.  Bui  Id. m 
calls  missile. m  to  get  plant  state  space  information.  Using  the  transfer 
funct  ion  G  from  Chapter  7  and  the  chosen  weighting  filters,  the  automated 
and  H  design  functions,  h2opt2.m  and  hinf3.m,  are  called  to  form  the 
controllers.  If  requested,  Bode  plots  of  the  filters  as  seen  in  Chapter  7 
are  also  provided. 

Analyze. m  calls  missile. ra  to  get  needed  transfer  functions  and  then 
produces  the  requested  plots.  The  plots,  which  can  be  seen  in  Chapter  6, 
include  y  time  response  and  singular  value  plots,  noise  attenuatuat ion 
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designer 

inputs 


Figure  B. 1  Computer  Codes  Used  In  Controller  Design  And  Evaluation 

singular  value  plots,  and  p  analysis  plots.  Mureal.m  is  called  to  produce 
the  p  analysis  plots.  Analyze3.m  is  similar  to  analyze. m  except  it  performs 
the  analysis  for  three  controllers  at  once  for  comparison. 

Simulate. m  calls  missile. m  to  get  missile  state  space  information.  It 
uses  Equation  (3.7)  to  calculate  the  missile  simulation  transfer  function 
and  then  performs  the  simulation  as  shown  in  Figure  3.1.  If  the  user 
doesn’t  input  simulation  initial  conditions,  defaults  will  be  used.  The 
0.001  second  time  step  simulations  could  take  up  to  four  hours  Vax  CPU  time 
to  run  so  they  were  run  in  batch  mode.  Simulate. m  also  has  options  to 
provide  many  plots  in  addition  to  missile  flyout  for  use  in  diagnosing 
errors  and  confirming  reasonable  simulations. 
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Since  the  navigational  command  block  plays  an  important  role  in 
resulting  missile  miss  distances,  guide. m  was  written  to  study  the  effects 
of  changes  in  navigational  constants  and  differences  in  execution  lag  time 
independent  of  the  missile  dynamics.  Guide. m  is  basically  simulate. m  with 
all  missile  dynamics  stripped  out.  A  missile  g  limiter  can  be  used  to  keep 
the  flyouts  realistic. 
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